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This article discusses quantum fluctuation properties of a crystal lattice, and in particular, 
phonon squeezed states. Squeezed states of phonons allow a reduction in the quantum fluctuations 
of the atomic displacements to below the zero-point quantum noise level of coherent phonon states. 
Here we discuss our studies of both continuous-wave and impulsive second-order Raman scattering 
mechanisms. The later approach was used to experimentally suppress (by one part in a million) 
fluctuations in phonons. We calculate the expectation values and fluctuations of both the atomic 
displacement and the lattice amplitude operators, as well as the effects of the phonon squeezed states 
on macroscopically measurable quantities, such as changes in the dielectric constant. These results 
are compared with recent experiments. Further information, including pre prints and animations, 
are available in jittp: / /www-personal. engin.umich.edu/~nori/squeezed. html 



I. INTRODUCTION 

Classical phonon optics [|l| has succeeded in producing many acoustic analogs of classical optics, such as phonon 
mirrors, phonon lenses, phonon filters, and even "phonon microscopes" that can generate acoustic pictures with 
a resolution comparable to that of visible light microscopy. Most phonon optics experiments use heat pulses or 
superconducting transducers to generate incoherent phonons, which propagate ballistically in the crystal. These 
ballistic incoherent phonons can then be manipulated by the above-mentioned devices, just like in geometric optics. 
Phonons can also be excited phase-coherently. For instance, coherent acoustic waves with frequencies of up to 10 10 
Hz can be generated by piezoelectric oscillators. Lasers have also been used to generate coherent acoustic and optical 
phonons through stimulated Brillouin and Raman scattering experiments. Furthermore, in recent years, it has been 
possible to track the phases of coherent optical phonons |2j due to the availability of femtosecond-pulse ultrafast 
lasers (with a pulse duration shorter than a phonon period) , and techniques that can measure optical reflectivity with 
accuracy of one part in 10 6 . In most situations involving phonons, a classical description is adequate. However, at low 
enough temperatures, quantum fluctuations become dominant. For example, a recent study [gj shows that quantum 
fluctuations in the atomic positions can indeed influence observable quantities (e.g., the Raman line shape) even 
when temperatures are not very low. With these facts in mind, and prompted by the many exciting developments in 
classical phonon optics, coherent phonon experiments, and (on the other hand) squeezed states of light we would 
like to explore phonon analogs of quantum optics. In particular, we study the dynamical and quantum fluctuation 
properties of the atomic displacements, in analogy with the modulation of quantum noise in light. Specifically, we 
study single-mode and two-mode phonon coherent and squeezed states, and then focus on a polariton-based approach 
to achieve smaller quantum noise than the zero-point fluctuations of the atomic lattice. The concepts of coherent 
and squeezed states were both originally proposed in the context of quantum optics. A coherent state is a phase- 
coherent sum of number states. In this state, the quantum fluctuations in any pair of conjugate variables are at the 
lower limit of the Heisenberg uncertainty principle. In other words, a coherent state is as "quiet" as the vacuum 
state. Squeezed states Q are interesting because they can have smaller quantum noise than the vacuum state in 
one of the conjugate variables, thus having a promising future in different applications ranging from gravitational 
wave detection to optical communications. In addition, squeezed states form an exciting group of states and can 
provide unique insight into quantum mechanical fluctuations. In recent years, squeezed states are also being explored 
in a variety of non-quantum-optics systems, including ion-motion and classical squeezing molecular vibrations 
||, polaritons @||, and phonons in crystals |p|-|l3|. Ref. 10-1^] propose a second-order Raman scattering (SORS) 



process for phonon squeezing: if the two incident light beams are in coherent states, the phonons generated by the 
SORS are in a two-mode squeezed state. Here we first present an introduction to the subject of coherent and squeezed 
phonons, and later on consider both the continuous wave case studied in [Io|Jl^| and the impulsive case studied in 
[12,E3. Squeezed phonons could be detected by measuring the intensity of the reflected or transmitted probe light 
[10 ]l2| , p^[ . This method has been used to detect coherent phonon amplitudes, since reflectivity and transmission are 
closely related to the atomic displacements in a crystal. Measuring a transmitted probe light pulse, Ref. [H observed 
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squeezed phonons produced by an impulsive SORS. The intensity of the CW SORS signal for many materials might 
be too weak to be detected with current techniques, but might be accessible in the future. 

II. ANALOGIES AND DIFFERENCES BETWEEN PHONONS AND PHOTONS 

Coherent and squeezed states were initially introduced to describe photons. Here we are interested in applying 
these concepts to phonons. Although both photons and phonons are bosons, they do have significant differences, and 
the physics of squeezed states of light cannot be straightforwardly extended to phonons. Table I is a chart presenting 
a schematic comparison between phonons and photons. Below we briefly mention a few important similarities and 
differences that are relevant to our study. 

Photons are elementary particles with no internal structure, thus are sometimes called simple bosons. On the other 
hand, phonons describe the collective displacements of very many atoms in a crystal, and are thus sometimes described 
as composite bosons fl5|| . Phonons are bosons because of the commutation relation between the coordinate and 
momentum operators. Kohn and Sherrington [fil| pioneered the research on composite bosons like phonons, excitons, 
etc., and classified them into two categories, with type-I referring to those bosons composed of an even number of 
fermions (such as 4 Hc atoms), and type-II referring to those that are collective excitations — such as phonons, excitons, 
magnons, etc. In this sense it is also possible to consider photons to be type-II composite bosons Jl5), because they are 
the energy quanta of electromagnetic field modes. Their commutation relation originates from the simple harmonic 
oscillators that are used to quantize the electromagnetic field. Essentially, both photons and phonons are field quanta: 
photons are quanta of a continuous field, while phonons are quanta of a discrete field. Non-interacting phonons are 
used to describe harmonic crystal potentials. However, anharmonicity, which leads to phonon-phonon interactions, 
is always present. Some properties of solids, such as lattice heat conductivity and thermal expansion, solely depend 
on the anharmonic terms in the crystal potential. In other words, phonons in general interact with each other. For 
photons, the situation is somewhat different. In vacuum and at low intensity, photon interactions are so weak that 
the rule of linear superposition holds. However, in nonlinear media, photons are effectively interactive, with their 
interaction mediated by the atoms. As mentioned above, phonons exist in discrete media. Therefore, phonons have 
cut-off frequencies, which put an upper-limit to their energy spectra. For a diatomic lattice, this limit is of the order of 
0.1 eV, which is in the infrared region. Photons, on the other hand, do not have such an upper bound for their energy. 
In addition, the discrete atomic lattice and the massive atoms lead to a finite zero-point fluctuation in the phonon 
field, while the continuous photon modes and the massless photons contribute to a divergent zero-point fluctuation 
in the photon field. The dispersion relations for photons and phonons are qualitatively different. Photons in free 
space have a linear dispersion relation. On the other hand, phonons have complicated nonlinear dispersion relations 
which generally have several acoustic and optical branches. The acoustic branches are linear around the center of the 
first Brillouin zone, i.e., the k = point, which is at the continuum limit. When the quasi-wave- vector k is close to 
the first Brillouin zone boundary, to saturates. The optical branches of the phonons have a different profile. Their 
dispersion relations are flat near k = 0, where to — loq. Furthermore, as k increases, u> decreases; indeed, the optical 
phonon dispersion relation can be even more complicated depending on the lattice structure. Compared to photons, 
which generally have relatively simple dispersion relations, phonons have nonlinear dispersion relations that make it 
more difficult to satisfy both energy and momentum (in fact, quasi- momentum) conservation laws simultaneously. 
The order of magnitude of the crystal cohesion energy determines that phonons have very low energies. In addition, 
phonons can easily couple to many other excitations which are in a similar energy range, and be perturbed by thermal 
fluctuations even at low temperatures. All these couplings make phonon dynamics very dissipative. Due to its strong 
damping, coherent phonons have very short lifetimes (~ 50 picosecond for optical phonons, while larger for acoustic 
phonons) H. On the other hand, there exist many materials in which photons can propagate with little dissipation, 
and furthermore very long photon coherent times can be produced by lasers. To summarize this brief comparison, 
we notice that the differences between phonons and photons can often hinder our effort to apply ideas originating in 
quantum optics to phonons. For example, there does not exist good phonon cavities at the moment, so that it is very 
difficult to select phonon modes. For phonons, we have to almost always deal with a continuum of phonon modes. In 
addition, short phonon lifetimes and the strongly dissipative environment of phonons also complicate the problems. 
These will be taken into consideration when we work on the theory. 
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III. PHONON OPERATORS AND THE PHONON VACUUM AND NUMBER STATES 



A phonon with quasimomentum p = Kq and branch subscript A has energy e qA = huj^x; the corresponding creation 
and annihilation operators satisfy the boson commutation relations 

Va'i b q\ = <W<5aa' , [bqx, Va'] = 0- (!) 
The atomic displacements Ui a of a crystal lattice are given by 

N 



qA 

Here refers to the equilibrium lattice positions, a to a particular direction, and <2 q is the phonon normal-mode 
operator 



^ = ^(M + ^a)- (3) 

For simplicity, hereafter we will drop the branch subscript A, assume that t/ qQ is real, and define a q-mode dimen- 
sionless lattice amplitude operator: 

u(±q) = 6 q + 6 f _ q + 6_ q + & q . (4) 

This operator contains essential information on the lattice dynamics, including quantum fluctuations. It is the phonon 
analog of the electric field in the photon case. When no phonon is excited, the crystal lattice is in the phonon vacuum 
state |0). The expectation values of the atomic displacement and the lattice amplitude are zero, but the fluctuations 
will be finite: 



((Au iQ ) 2 ) vac = ((u M ) 2 )vac ~ (^a)vac ( 5 ) 



/vac 

N 



^ 2NmLu qa ' V ' 

q 

((Au(±q)) 2 ) vac -2. (7) 

The eigenstates of the harmonic phonon Hamiltonian are number states which satisfy 6 q |n q ) = ^/n^n^ — 1). The 
phonon number and the phase of atomic vibrations are conjugate variables. Thus, due to the uncertainty principle, the 
phase is arbitrary when the phonon number is certain, as it is the case with any number state |n q ). Thus, in a number 
state, the expectation values of the atomic displacement (n q |ztj Q |n q ) and q-mode lattice amplitude (n q |u(±q)|n q ) 
vanish due to the randomness in the phase of the atomic displacements. The fluctuations in a number state |n q ) are 

((Au ia )\ um = + jr , (8) 



q'/q 



((Au(±q)) 2 ) num = 2 + 2n q . (9) 



IV. PHONON COHERENT STATES 

A single-mode (q) phonon coherent state is an eigenstate of a phonon annihilation operator: 

& q |/3 q ) = /3 q |/3 q ) . (10) 
It can also be generated by applying a phonon displacement operator _D q (/3 q ) to the phonon vacuum state 

|/3 q ) = £>q(/? q )|0) - exp(/3 q 6 q - /3 q 6 q )|0) (11) 

exp ( 9 ) 1^ "7^7 |n <^ • (12) 



n a =o v n q- 
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As is shown above, a phonon coherent state is a phase-coherent superposition of number states. Moreover, coherent 
states are a set of minimum-uncertainty states which are as noiseless as the vacuum state. Coherent states are also 
the set of quantum states that best describe the classical harmonic oscillators |t6). A single- mode phonon coherent 
state can be generated by the Hamiltonian 

H = hco^ (b^ + ~) + A q (i) 6 q + A q (i) 6 q (13) 

and an appropriate initial state. Here A q (<) represents the interaction strength between phonons and the external 
source. More specifically, if the initial state is a vacuum state, | -0(0)} = |0), then the state vector becomes a single-mode 
coherent state thereafter 

|^)) = |A q (i)e-^), (14) 

where 

A,(*) = -ijT Vrje^dr (15) 

is the coherent amplitude of mode q. If the initial state is a single-mode coherent state |^(0)) = |a q ), then the state 
vector at time t takes the form 

|V>(i)> = |{A q (i)+a q }e-^ (16) 

which is still a coherent state. In a single-mode (q) coherent state |A q (t) e~ luJcit ), (ui a (t))coh and (u(±q)) co h are 
sinusoidal functions of time. The fluctuation in the atomic displacements is 

q 4 

The unexcited modes are in the vacuum state and thus all contribute to the noise in the form of zero point fluctuations. 
Furthermore, 

((Au(±q)) 2 ) coh = 2. (18) 

^From the expressions of the noise ((Aui Q ) 2 ) co h and ((Au[±q)] 2 ) co h, it is impossible to know which state (if any) has 
been excited, while this information is clearly present in the expression of the expectation value of the lattice amplitude 
(zt(±q)) co h- These results can be straightforwardly generalized to multi-mode coherent states. Coherent phonons have 
been the subject of considerable interest in recent years. Typically, the dynamics of coherent phonons are described 
by using classical equations of motion. Here we present a quantum description and show that it is consistent with the 
classical one and, as an additional bonus, contains information on quantum fluctuations. Coherent phonons can be 
generated by a femtosecond short pulse laser [ p!o[ . A femtosecond pulse duration is much shorter than any phonon 
period and therefore acts as a delta-function driving force. It can produce coherent longitudinal optical (LO) phonons 
P,[To|. We can make a very simplified calculation by replacing the coupling strength A q (i) with A6(t — to) in Eq. (|l3|), 
so that 

H coh = fiw q 6t ft q + A 6(t - t ) 6 q + A* 5(t - t ) 6 q . (19) 

Here A = \A\e z ^ A is a time- independent complex amplitude containing the information of the photon-phonon inter- 
action and the coherent amplitude of the relevant modes in the incident optical pulse. We assume that the crystal is 
in the phonon vacuum state before it is hit by the laser pulse at t — to- The time-evolution operator U(t,to) for the 
phonon mode can be written as 

U(t, to) = exp (-iw q t 6 q 6 q ) exp f-i^e-*** - ^& q e^*°) • (20) 



In other words, for t > to, the crystal is in a coherent state |A q e The coherent phonon amplitude A q is 

iA* 
~2h 



A q = -^e^ t0 . (21) 
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which is a constant complex number. Thus, a very short laser pulse conveniently provides a time-independent 
amplitude and a coherent phase to the active phonon mode(s). For t > to, the average of the atomic displacement 
operator in the state jAqe -4 ""') ® |0_ q ) becomes 



2h 



Nmujq 



|[/ qQ A q | mx[(u q (t - t ) - q • R; - <j> A - <j>u) ■ (22) 



These longitudinal optical phonons can have a coherence time of about 50ps at 10K, and even longer at lower 
temperatures In the classical sense, "coherent" means a wave with a well-defined phase, or waves that can 
interfere with each other when superimposed. Here we have shown that a single-mode coherent state of phonons 
generated by a short laser pulse is indeed a plane wave with a well-defined phase. Thus, these phonons in a quantum 
coherent state are also coherent in a classical manner. Furthermore, they are in a minimum-uncertainty state. 



V. PHONON SQUEEZED STATES 

In order to reduce quantum noise to a level below the zero-point fluctuation level, we need to consider phonon 
squeezed states. Quadrature squeezed states are generalized coherent states |p7| . Here "quadrature" refers to the 
dimensionless coordinate and momentum. Compared to coherent states, squeezed ones can achieve smaller variances 
for one of the quadratures during certain time intervals and are therefore helpful for decreasing quantum noise. 
Figures [l] and || schematically illustrate several types of phonon states, including vacuum, number, coherent, and 
squeezed states. These figures are the phonon analogs of the illuminating schematic diagrams used for photons Jl7| . 
A single-mode quadrature phonon squeezed state is generated from a vacuum state as 

|a q ,0 = I>q(aq)5 q (0|0); (23) 
a two-mode quadrature phonon squeezed state is generated as 

l^qi ! Q; q2 ' 

= A,> q JA,> q2 )S qi , q2 (£)|0) . (24) 
Here _D q (a q ) is the coherent state displacement operator with a q = |a q |e i<6 , 



S q (0 = exp ( t-b\ - |fe q 2 ) , (25) 



S qi , q2 (0 = exp - &\b\ 2 ) , (26) 

are the single- and two- mode squeezing operator, and £ = re t9 is the complex squeezing factor with r > and 
< < 2t:. The squeezing operator 5' qiiq2 (5) can be produced by the following Hamiltonian: 

-frqi,^ = ftUq^b^ + 7kj q2 6 q2 6 q2 

K(^ qi & q2 +CWM q2 - (27) 
The time-evolution operator for such a Hamiltonian has the form 

U(t) = exp (~H ?j exp [f (t)6 qi 6 q2 - > ( 28 ) 

where 

H Q = nw qi & qi & qi + 7iw q2 6 q2 6 q2 , (29) 



= ((rj^Ki+^l^T. (30) 



— oo 



Here £(£) is the squeezing factor and £(i) is the strength of the interaction between the phonon system and the 
external source; this interaction allows the generation and absorption of two phonons at a time. The two-mode 
phonon quadrature operators have the form 
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X(q, -q) = 2- 3 / 2 (b q + &t + fe_ q + jQ (31) 
= 2- 3 / 2 U (±q), 

P(q, -q) = -*2~ 3/2 (t q - & q + 6_ q - fe f _ q ) . (32) 

We have considered two cases where squeezed states were involved in modes ±q. In the first case, the system is in a 
two-mode (±q) squeezed state |a q , a_ q , £), (£ = re 10 ), and its fluctuation is 

<[A W (±q)] 2 ) = 2 (V 2 ^ cos 2 ? + e 2r sin 2 . (33) 

In the second case, the system is in a single-mode squeezed state |a q ,£) (a q = |a q |e l<ii ) in the first mode and an 
arbitrary coherent state |/3- q ) in the second mode. The fluctuation is now 

([Au(±q)] 2 ) = l + e 2r sin 2 



-e 



-2r „ 2 



2 

In both of these cases, ([Au(±q)] 2 ) can be smaller than in coherent states (see Fig. 2) 



cos 2 U+- . (34) 



VI. SECOND ORDER RAMAN SCATTERING (SORS) 

So far we have focused on concepts and the actual quantum mechanical states relevant to coherent and squeezed 
phonons. Now we will focus on one particular approach to generate squeezed phonon states via second order Raman 
scattering. The SORS process originates from the quadratic term in the polarizability change 5P a p of a crystal. The 

photon-phonon interaction V that leads to the SORS process is Q V — — \ J2 a p E q £«' ^ad ^ QtdQ-qj'EiaEzp ■ 
Here, E\ a and E^p are electric field amplitudes along a and (3 directions with frequencies u>± and u>2- The second- 
order polarizability tensor -P^g _qL5 satisfies P^a^ = P^ 3 ™ ~ ^aft*^ • -R- ecan that the complex phonon normal 
mode operator Qqj of the phonons is related to the phonon creation b'L^ and annihilation bqj operators by Qqj = 

bqj + frLqj- If the incident photon fields are not attenuated we can treat the optical fields as classical waves, and also 
consider the different pairs of ±q modes as independent, and treat them separately. Thus, for one particular pair 
of ±q modes, the complete Hamiltonian for the two phonon modes involved in the SORS process has the form |L8|| : 
Hq = Hq — {4 _1 J2af3 P^jT^ElaEspyQqQ-q , where Hq — hujqlb^bq + b^ q &- q } is the free phonon Hamiltonian for 
the modes q and — q, aj q = (uj\ — W2)/2, and the branch labels j and j' have been dropped. Here we consider two 



different cases. The first is when the incident photons are in two monochromatic beams [|10 12 1; i.e., with electric fields 
Ej = £j cos(ujjt + <pj); j = 1, 2. In the second case the incident photons are in an ultrashort pulse whose duration is 
much shorter than the phonon period |l^ , |i"4| ]. 



VII. SQUEEZED PHONONS VIA CONTINUOUS WAVE SORS 

Let us now first consider the continuous wave (CW) case. Because the photons are monochromatic, we can take a 
rotating wave approximation pSj and keep only the on-resonance terms in the Hamiltonian. The off-resonance terms 
only contribute to virtual processes pQ] at higher orders. This approximation is appropriate for times much longer 
than the phonon period. The simplified Hamiltonian has the form 

Hq CW) =Hq-Xq { bqb^q e 2 -V+#12 + ) 

A q = ^iE p 5r q " 2 ^ ( 35 ) 

a/3 

where (^12 and A q refer to the overall phase and amplitude, respectively, of the product of the 2nd-order polariz- 
ability and the incident electric fields. Recall that P^p~ cl is real, therefore the phase 0i2 has no q-dependence. It 
originates solely from the two photon modes. The Schrodinger equation for the ±q-mode phonons is ifidt \ipq{t)) = 



6 



TCq (t) \?pq(t)) , and its time-evolution operator can be solved by a transformation into the interaction picture. The 
result can be expressed as |l0| , [l^] 

|^ q (t)) = e^t/ih} e {c;M- q -C q b^L q } |^ q ( 0)> ; ( 36) 

where £ q = — iXq t e~ 1( ^ 12 /h . Notice that the second factor in the time-evolution operator is a two- mode quadrature 
squeezing operator p]j |. In the CW case considered here, the amplitude of the squeezing factor £ q grows linearly 
with time. However, this initial linear growth will be eventually curbed by subsequent phonon-phonon scattering 
and optical pump depletion. In other words, the expression for the squeezing factor £ q is valid for times much larger 
than one phonon period, but much smaller than phonon lifetimes (because this treatment considers non-decaying 
phonons). By solving the phonon Langevin equation, we have shown that the squeezing factor in the CW SORS will 
eventually saturate at a constant value determined by the strength of SORS and the phonon decay constant In 
addition, the phase of the squeezing factor is determined by the phase difference of the two incoming light waves. If 
the ±q phonon modes are initially in a vacuum state or in a coherent state, the SORS will drive them into a two-mode 
quadrature squeezed state fl0| , |l2|| . The time evolution operator of all the phonon mode pairs (instead of just one 
pair of ±q modes) that are involved in this SORS process has the form U(t) = JT t/ q (t) . Therefore, as long as the 
photon depletion is negligible, all the phonon modes that are involved in a SORS process are driven into two-mode 
quadrature squeezed states. In other words, squeezing can be achieved in a continuum of phonon modes by a CW 
stimulated SORS process. 

VIII. SQUEEZED PHONONS VIA IMPULSIVE SORS 

Recently, an impulsive SORS process has been used to experimentally generate phonon squeezing |Q . Here we treat 
the problem expressing the time evolution operator of the system in terms of a product of the two-mode quadrature 
squeezing operator and the free rotation operators |23|| . Since the incident photons are now in an ultrashort pulse, the 
complete Hamiltonian can be solved in the limit when the optical field can be represented by a <5-function. Such an 
approximation is usually considered when the optical pulse duration is much shorter than the optical phonon period, 
which is experimentally feasible with femtosecond laser pulses. The Hamiltonian for the SORS can now be written as 
Tt = J2q {-^q — -^WQqQ-q} ; where A q carries the information on the amplitudes of the incoming optical fields 
and the electronic polarizability. Notice that the light-phonon coupling strength A q in the CW case has units of 
energy, while A q here has units of h. To further simplify the problem, we assume that only ±q modes are involved in 
the process. Such a simplification is possible when the photon depletion and the phonon anharmonic interaction are 
negligible, so that different pairs of phonon modes are independent from each other. The Hamiltonian is now 

K = ff q -A q ^)Q q Q_ q , (37) 

and the Schrodinger equation for these two phonon modes is ihdt\ipq{t)) = 7i q \ipq(t)}. This equation can be solved by 
separating the free oscillator terms and the two-phonon creation and annihilation terms. The resulting time-dependent 
wavef unction is 

X exp {C q *M- q - C'qbqblq}\lpq(0-)) • (38) 

Here £ q = — iA q e^ lX ^ h jh. Hence the effect of the optical pulse is clear: it first applies a two-mode quadrature 
squeezing operator on the initial state, then rotates the state by changing its phase p3]| . The state will then freely 
evolve after t = + . This result is consistent with Ref . fl4|| where the time-evolution operator is expressed in terms of 
real phonon normal mode operators ^M , instead of the complex ones used in this paper. Notice that, in contrast to 
the CW SORS, the phase of the squeezing factor (' for the impulsive case is fixed by the intensity of the light pulse. 

IX. MACROSCOPIC IMPLICATIONS AND TIME-DEPENDENCE OF THE DIELECTRIC CONSTANT 

Now that we have obtained the phonon states for both the CW and pulsed SORS cases, let us consider the 
macroscopic implications of these states. More generally, let us first discuss the implications of the phonon squeezed 
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states disregarding how they are generated. An experimentally observable quantity O which is related to the atomic 
displacements in the crystal can generally be expressed in terms of Q q : O = O(0) + J2qi^O/dQq)Qq + . . . = Oo + 
0\ + 2 + ■ ■ ■ where the first term Oo = 0(0) is the operator O when all Q q 's vanish. An example of an experimentally 
observable quantity O is the change in the crystal dielectric constant Se due to the atomic displacement produced by 
the incident electric fields. To first order in Q q , Se = Sex = E 9 ,>o + 6 -q) e ** q + ( b ~i + 6 q) e ~ 4 * q ] ' 

Here * q is the phase of dO/dQq — d(Se)/dQq. Indeed, a widely used method to track the phases of coherent 
phonons in the time-domain |2| is based on the observation of the reflectivity (or transmission) modulation SR (ST) 
of the sample, which is linearly related to Se — the change in the dielectric constant due to lattice vibrations. The 
above equation for Se indicates that we can define a generalized jn| lattice amplitude operator PJH[|: u s (±q) = 

(bq + foLq^ e l *i + (6_ q + b q ) 6^**1 . This generalized lattice amplitude u s (±q) = 2i?e{Q q e l *i} is the underlying 
microscopic quantity related to an observed reflectivity or transmission modulation when the linear term in Q q , Se±, 
exists. Even if (Se\) vanishes, (A(Sei) 2 } does not. Since different pairs of ±q phonon modes are uncorrelated to one 
another, the fluctuation of Se\ can be expressed as ((A<5ei) 2 ) = E fe >o ( h / 2 ^) \d(Se) /9Q q | (Au|(±q)) . Here the 
state is \ip(t)) — U(t)\ip(0)) = n q ^q(*)l^q(0)) m either the CW or the impulsive case. We can again focus on a single 
pair of ±q modes. In the CW case, using Eq. (|36|), the fluctuation is 

(Au 2 g (±q))^ = 2{e^cos 2 (fi q (t) + 12 /2) 

+e 2r "sin 2 (a i (i) + 12 /2)}, (39) 

where r q = |£ q | = A q t/h, flq(t) — ojqt + 7r/4, and hereafter (. . .) denotes an expectation value on squeezed states, 
unless stated otherwise. Therefore, at certain times, the fluctuation (Au 2 (±q))( cl ") can be smaller than 2, which is 
the vacuum fluctuation level. Furthermore, all the pairs of phonon modes that are driven by the stimulated SORS 
process share the same frequency: ujq — (uji — w 2 )/2. Therefore, all the fluctuations (Aw 2 (±q))( c "') evolve with the 
same ujq. Notice that there is no dependence on ^ q in the final expression of (Au 2 g (±q))( cw \ and the squeezing factor- 
phase 012/2 has no q-dependence, all the pairs of modes involved through the SORS share the same phase in their 
fluctuations. Therefore there can be squeezing in the overall fluctuation ((ASei) 2 )^ cw K Furthermore, the phase of this 
overall fluctuation can be adjusted by tuning the phase difference of the two incoming light beams. In the impulsive 
case if the ±q-mode phonons are driven into a squeezed vacuum state, the fluctuation in u g (±q) is 

(Aw 2 (±q))' = 2{e~ 2 < cos 2 %{t) + e 2 < sin 2 O q (i)} , (40) 

where r q = = A q /?i, and fi q (£) = Qq(t) — r q . Again, the squeezing will reveal itself through oscillations in 
([A((5ei)] 2 (q))' which is proportional to (Au 2 (±q))'. Note that these oscillations are essentially the same as the ones 
obtained in the CW case. However, now the squeezing factor is time-independent. Also, the t — phase 7r/4 — r q in 
Eq. j40| ) is q-dependent. Eq.(pO[) can be rewritten as (Au 2 (±q))' = 2 {cosh2r q + sinh2r q sin(2oj q t — r q )} . For small 
r q , this becomes (Ait 2 (±q))' = 2{1 + 2r q 2 + 2r' sin(2cj q i — f q )}- This expression has essentially the same form as the 
one obtained in @: (Q 2 (£)) = (O q (0)){l + 2^ + 2£ q sin(2w q £ + ip q )}. The small phase term ipq is neglected in Q 
when computing transmission changes. The difference in phases, r q versus </? q , is negligible in the limit of very small 
squeezing factor, and originates from the different interaction Hamiltonians used here and in [pL The interaction term 
i n fli] ] is proportional to u 2 (±q) with VPq = (notice that their Q q is real and based on standing wave quantization 
Jl8|). Therefore, the interaction Hamiltonian in [|l4| is (in our notation) V cx w 2 (±q) cx 20 q O_ q +0 2 +Q 2 _ q . However, 
the last two terms in this expression do not satisfy momentum conservation, we thus did not include them and kept 
only QqQ-q in our interaction term (this form is also used by Ref. [^8|). When the linear perturbation Sei due to 
phonons vanishes, such as in fl4|] , the second order correction 2 (= <5e 2 ) must be considered. When the phonon states 
are modulated by a SORS, so that the ±q modes are the only ones which are correlated, Se 2 — J2q dQ dQ — QqQ-q ■ 
Let us first focus on one pair of ±q modes in the CW case. In a vacuum state, (0|O q O- q |0) = 1 ; while in a squeezed 
vacuum state |0) sq , sq (0|Q q O- q |0) sq = (Au 2 (±q))(™)/2, with the right hand side given in Eq. (|39"[). Therefore, the 
expectation value of Q q 0- q in a squeezed vacuum state is periodically smaller than its vacuum state value. Let us 
now include all the phonon modes that contribute to (5e 2 . In a vacuum state, (0|5e 2 |0) = ^ q d 2 Se/ (dQq 9Q_ q ) . On 
the other hand, in a squeezed vacuum state, 

w-sEs^**^'- (41) 

Since the phase </>i 2 /2 has no q-dependence, contributions from the phonon modes sharing the same frequency add 
up constructively. It is thus possible that (5e 2 ) is periodically smaller than its vacuum state value. Similarly, in 
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the impulsive case, (Se^)' = 2 1 X) q dQ 3Q — (Aftj(±q))'; However, the phase factor in {Se^) 1 has a q-dependence 
through r q , so that all the phonon modes with the same w q do nof contribute to (662)' synchronously. In the CW 
SORS and in the very-small-r q limit impulsive SORS the phase of the expectation value (QqQ-q) does not depend on 
q; this is crucial to the experimental observation of modulations in the dielectric constant, because this q-insensitivity 
leads to constructive summations of all the q pairs involved. Also, at a van Hove singularity a large number of modes 
contribute to 5t2 with the same frequency and phase, thus their effect is larger and easier to observe 14 1. 



X. SQUEEZED PHONONS VIA A FINITE- WIDTH SORS 



Real light pulses are not (^-functions. Therefore, we have also considered a SORS pumped by a light pulse with a 
finite width (smaller than the phonon period T) instead of a (^-function. For a fixed peak height /, we find |22| that 
the optimal pulse width T° pt that maximizes the squeezing effect satisfies T° pt w T/4.4. This calculation indicates 
that the experiments (TJ] used a pulse width which is nearby the optimal value (T/4.4 w 300/4.4 fs « 68 fs ~ T p ). 
The calculation can be summarized as follows. First, in the impulsive Hamiltonian we replace the (5-function 
by a Gaussian with its width T p as a variational parameter. Since now the Hamiltonian is time-dependent in the 
interaction picture, we cannot directly integrate the Schrodingcr equation. Instead, we use the Magnus method to 
obtain the time evolution operator and keep only the dominant first term. This approximation is valid when the pulse 
duration is shorter than the phonon period. We then calculate the width T° pt of the Gaussian that maximizes the 
squeezing factor. For a constant peak intensity, a pulse that is too narrow does not contain enough photons; while it 
can be proven that a pulse which is too long (i.e., with a width comparable to T), attenuates the squeezing effect. 



XI. PHONON SQUEEZING MECHANISM 



What is the mechanism of phonon squeezing in the SORS processes? For the CW case, the Hamiltonian is the same 
as an optical two-mode parametric process Jill , with the low frequency interference of the combined photon modes as 
the pump, the two phonon modes as the signal and idler. The frequencies of these modes satisfy cj q + cj_ q = uji — oj 2 - 
The impulsive case is slightly different. Although the Hamiltonian is similar to a parametric process, the energy 
transfer from the photons to the two phonon modes is instantaneous. The resulting phonon state is a two-mode 
quadrature squeezed vacuum state. Indeed, a regular parametric process pumps energy into the signal and idler modes 
gradually, while the impulsive SORS does it suddenly. The correlation between the two phonon modes, and thus the 
squeezing effect, is also introduced instantaneously Notice that this mechanism is reminiscent of the frequency-jump 
mechanism proposed in B. In the impulsive SORS, the frequency of the phonon modes has an "infinite" <5-peak 
change at t = 0, while the 

frequency-jump mechanism has finite frequency changes, and squeezing there can be intensified by repeated fre- 
quency jumps at appropriate times. However, as it has been pointed out in Q, a finite frequency jump up immediately 
followed by an equal jump down results in no squeezing at all. 



XII. EXPERIMENTAL SEARCH FOR SQUEEZED PHONONS 



One group [Q has claimed to detect squeezed phonons. Regrettably, no other group has succeeded in reproducing 
their results. Recent calculations indicate that these initial experiments [fl4| do not prove that phonon squeezing 
has been achieved experimentally. Here we discuss several important points to be considered for the eventual future 
experimental observation of a phonon squeezed state. First of all, squeezing must refer to a phonon mode with a 
variance that falls below the standard quantum limit. This crucial comparison has not been produced yet. Thus, a 
reliable determination of the vacuum noise level is necessary in order to establish if a squeezing of the vacuum noise 
has indeed occurred. Experiments in quantum optics have an independent way to reliably obtain the noise level of the 
vacuum: by using an independent light beam as a local oscillator. Currently, there is no such phase-sensitive detection 
scheme for phonons. However, as we point out later, a pump-probe experiment may be able to establish such a criteria, 
but discretion has to be applied. Second, for phonons, with relatively low energy compared to photons, thermal noise 
should always be considered in an experiment. This is especially important when the noise modulation factor is small, 
such as the case in the reported experiment |L4]], where the noise modulation factor is only 0.0001%. To achieve total 
noise modulation, as can be done through CW and impulsive SORS, is not equivalent to the squeezing of the vacuum 
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noise. Only when the noise modulation factor is big enough to overcome the thermal noise, then the quantum noise is 
suppressed. Third, in that experiment, the squeezing factor is obtained from the modulation of the light transmission 
through a crystal. However, this change in light transmission describes the modulation of the total fluctuations (both 
quantum and thermal) of the atomic displacements. Thus, by itself this modulation of the noise is not a proof of 
quantum noise suppression, let alone of squeezing below the vacuum noise level. Fourth, the efficiency of the signal 
detection is a crucial issue in photon squeezing experiments, but not addressed in its phonon "analog" experiment. 
In a homodync or heterodyne detector for photons, if the efficiency of the photon counters is less than one, additional 
noise is introduced into the signal. Similarly, in a pump-probe phonon detection scheme, the probe light pulse and 
photodetector should introduce additional noise into the final signal. Therefore, a careful analysis of these additional 
noise sources should be performed when employing a pump-probe scheme to detect phonon squeezing. For the above 
reasons, it is premature and unwarranted to claim that squeezed phonons have been observed experimentally. The 
evidence presented so far is incomplete and inconclusive. The answer to the question: "Can phonons be squeezed like 
photons?" is yes on theoretical grounds, but the experimental proof still lies in the future. The rest of this section 
presents quantitative derivations in support of the statements made above. Even with the reservations presented 
above regarding the pump-probe scheme, we would like to point out that it does provide a possible means for the 
detection of phonon squeezing. Here we discuss the criteria for achieving squeezing of quantum noise through an 
impulsive SORS process in a pump-probe experiment. We then give a numerical estimate on whether squeezing of 
quantum noise has been achived in Experimentally, the observed quantity is the change in the transmission T 
due to the impulsive SORS process. Up to second order in <3 q , T can be expressed as 

8T B 2 T 

where the average is over the phonon states of the crystal. In a squeezed vacuum (or thermal) state considered 
here here, (Q q ) = and (Q q Q_ q ) = 2(Au 2 (±q, <)). Therefore, T does contain informations on the atomic lattice 
fluctuations. Through impulsive SORS, we can achieve modulation of the total noise of the system, and the effect 
can be seen through a modulation of the transmission T, or more specifically, AT/T. Starting from a thermal state, 
we have: 

AT 1 ^ d 2 T 



f E nnnn [1 + 2iV K)K sin(2c q t) . 



T dQ^OQ 

Notice that this modulation of the total noise is not equivalent to a squeezing of quantum noise, because AT = is 
an indication of total noise limit (That is, quantum noise plus thermal noise), not quantum noise limit. Even though 
the experiment is done at low temperatures, so that thermal noise is miniscule, we still have to compare the thermal 
noise to the squeezing factor (or modulation factor through impulsive SORS) because the later is very small, too. 
One criteria for quantum noise suppression would be 

AT' 



T 



q 



q iy^_q 



A negative AT' /T here indicates that, at least in some spectral regions, the atomic displacement noise is squeezed 
below the ground state limit. We are now comparing the thermal population 2N(uj <i ) and the squeezing factor r q . In 
summary, there is always modulation of total noise by the impulsive SORS, but only when r q > 27V (cj q ) do we achieve 
noise smaller than the ground state limit. For example, at an experimental temperature of 10K, the corresponding 
thermal energy ksT is about 0.88meV. The acoustic phonons involved in the experiment reported in ||l4| have a 
frequency of about 2.7THz, which corresponds to an energy quantum of Tiw w 10.6meV. Therefore, the thermal noise 
factor here is about N = 1 / {exp(hui / hsT) — 1} w e~ 12 w 0.6 * 10 -5 , which is almost identical to the squeezing factor 
given in fl4f| . Therefore, it is quite clear that this particular experiment has not achieved squeezing of quantum noise 
below the vacuum noise limit. However, if the experimental temperature is further lowered, the thermal factor will 
become smaller. The squeezing action initiated by the second order Raman scattering should then be strong enough 
to realize suppression of noise below the ground state level. 
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XIII. CONCLUSIONS 



We have presented an overview of the definitions relevant to quantum phonon optics, including quantum coherent 
and squeezed phonons. Afterwards, we have studied theoretically the generation of phonon squeezing using a stimu- 
lated SORS process. In particular, we calculated the time evolution operators of the phonons in two different cases: 
when the incident photons are in monochromatic continuous waves, and when they are in an ultrashort pulse. The 
amplitude of the squeezing factor initially increases with time and then saturates in the CW SORS case, while it 
remains constant in the pulsed SORS case. In addition, the t = phase of the squeezing factor in the CW SORS, 
012, can be continuously adjusted by tuning the relative phase of the two incoming monochromatic photon beams, 
while for the pulsed SORS the phase (oc Aq) of the squeezing factor is determined by the amplitude of the incoming 
light pulse. For both cases we calculated the quantum fluctuations of a generalized lattice amplitude operator and the 
second order contribution to the change in dielectric constant, which is measurable. For the finite-width impulsive 
case, we computed the optimal pulse width, in terms of the phonon period, that maximizes the squeezing effect. 
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FIG. 1. Schematic diagram of the uncertainty areas (shaded) in the generalized coordinate and momentum (X(q, — q), 
■PCQi ~~ q)) phase space of (a) the phonon vacuum state, (b) a phonon number state, (c) a phonon coherent state, and (d) a 
phonon squeezed state. Here X(q, — q) and P(q, — q) are the two-mode (±q) coordinate and momentum operators, which are 
the direct generalizations of their corresponding single-mode operators. Notice that the phonon coherent state has the same 
uncertainty area as the vacuum state, and that both areas are circular, while the squeezed state has an elliptical uncertainty 
area. Therefore, in the direction parallel to the 9/2 line, the squeezed state has a smaller noise than both the vacuum and 
coherent states. 

FIG. 2. Schematic diagram of the time evolution of the expectation value and the fluctuation of the lattice ampli tude operator 
w(±q) in different states. Dashed lines represent (w(±q)), while the solid lines represent the envelopes (u(±q))±-^/ {[Aw(±q)] 2 ). 
(a) The phonon vacuum state |0), where (w(±q)) = and ([Au(±q)] 2 ) = 2. (b) A phonon number state |n q ,n_ q ), where 
(w(±q)) = and ([A«(±q)] 2 ) = 2(n q + n_ q ) + 2. (c) A single-mode phonon coherent state |a q ), where (u(±q)} = 2|a q | cosu q f 
(i.e., a q is real), and ((Au[±q)] 2 ) = 2. (d) A single-mode phonon squeezed state |a q e _l " q ', £(i)), with the squeezing factor 
£(t) = re' 2 ^^ and r = 1. Here, (u(±q)) = 2|a q | cosu; q i, and ([Au(±q)] 2 ) = 2(e~ 2r cos 2 w q i + e 2r sin 2 w q t). (e) A single-mode 
phonon squeezed state, as in (d); now the expectation value of u is («(±q)) = 2|a q | sincj q t, (i.e. a q is purely imaginary), and 
the fluctuation ([Au(±q)] 2 ) has the same time-dependence as in (d). Notice that the squeezing effect now appears at the times 
when the lattice amplitude (u(±q)) reaches its maxima, while in (d) the squeezing effect is present at the times when (u(±q)) 
is close to zero. 
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